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ABSTRACT 

Lyman Alpha Emitters (LAEs) are galaxies that have been selected on the basis of a 
strong Lya emission line in their spectra. Observational campaigns over the last decade 
have dramatically increased the sample of known LAEs, which now extends out to z = 
7. These discoveries have motivated numerous theoretical studies on the subject, which 
usually define LAEs in their models based on sharp Lya luminosity and equivalent 
width (EW) cuts. While broadly representative, this procedure does not mimic the 
selection from observational programs in detail, which instead use cuts in various 
colour-spaces. In this paper we investigate what implications this disjoint may have for 
theoretical studies that aim to model the observed population of LAEs. We construct 
an empirical model for the number density of star forming galaxies as a function of their 
UV and Lya luminosity, utilising measured constraints on the luminosity functions of 
drop-out galaxies, and their luminosity dependent probability distribution function of 
Lya EW. In particular, we investigate whether the LAE luminosity functions can be 
reproduced by defining LAEs using a (z-dependent) Lya luminosity and EW threshold. 
While we are able to reproduce the observed distribution of Lya EW among LAEs 
out to restframe EW~ 200 A, we find that our formalism over-predicts both the UV 
and Lya luminosity functions of LAEs by a factor of 2-3, and is inconsistent with 
observations at the ~ 95% level. This tension is partially resolved if we assume the 
Lya EW-distribution of drop-out galaxies to be truncated at restframe EW >150 

A. However the overprediction indicates that modeling LAEs with simple REW and 
luminosity cuts does not accurately mimic observed selection criteria, and can therefore 
lead to uncertainties in the predicted number density of LAEs. On the other hand, 
the predicted redshift evolution is not affected. We apply our formalism to drop-out 
galaxies at z > 6, and predict the luminosity functions of LAEs at z = 7 — 9. 

Key words: galaxies: high redshift - galaxies: high-redshift - galaxies: luminosity 
functions - line: formation - radiative transfer - scattering 



1 INTRODUCTION 

Two complementary observational techniques have been 
very succesful in finding high-redshift galaxies. The Lyman 
Break technique - or drop-out technique - has been used to 
constrain the observed rest-frame UV luminosity functions 
of Lyman Break galaxies (LBGs) out to redshifts a s high 
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iFinkelstein et al 



Bouwens et"afl 12010 j; iBunker et all l20ld: 



20ld : lYan et~ai1 l20ld ; lOesch et al.ll201ll 



The Lyman Break technique relies on the fact that the spec- 
tra of star forming galaxies are strongly suppressed blue- 
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ward of the Lyman-a resonance at A <J1216 A (|Steidel et al.l 
1996). The presence of this 'break' in the spectrum is a di- 
rect consequence of the temperature of stellar atmospheres, 
and of the interstellar and intergalactic absorption of ioniz- 
ing and Lyman series photons. The resulting break in the 
spectrum strongly suppresses the broad-band flux blueward 
of some-redshift dependent-filter that is observed from a 
galaxy (i.e. V drop-out galaxies at z ~ 5 are detected only 
in filters redder than V). 

Narrowband surveys have constrained luminosity func- 
tions of Lya e mitting galaxies ( a.k.a Lya emitt e rs, or LAEs) 



out to z ~ 9 llHu et all 12002 : live et all l200d; l Ouchi et al 



2008, l20ld : IWillis et al.H200Sl : 



Hibon et al.ll201d : iTilvi et al 
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120101 : IClement et al.ll2()ITh . while Lehnert et al. (2010) have 
also reported discovery of a Y 105 drop-out galaxy with strong 
Lya emission. The narrow band technique relies on the pres- 
ence of a strong Lya e mission line, which mostly originates 
in galactic HII regions (jPartridge fc PeeblesHl967r ). The pres- 
ence of such a line can produce an excess of observed flux in 
a narrowband filter (FWHM~ 100 A) that is larger than ex- 
pected bas ed on the observed fl ux in overlapping broadband 
filters (e.g. iRhoads et al.ll2000l h 

Existing observations have given us accurate determi- 
nations of the luminosity functions of LBGs and LAEs. 
The observed rest-frame UV luminosity function of drop- 
out gal axies decreases monotonically with redshift a t z ^ 
3 (e g. iReddv fc Steidell 120091 ; iBouwens et all 120061 . 120071 . 
I2008T ). These luminosity functions can be converted into 
a cosmic star formation rate density, which drops by 
more than an order of magnitude betwe en z = 3 and 



> 



6 (e.g. iHopkins fc Beacoml 120061 ; | Bouwens et alj 



120091 : iRobertson et al.ll2010t ). Furthermore, the broad band 



colours of drop-out galaxies become bluer with redshift, 
which indicates that star forming galaxies become increas- 
ingly dust-free dStanwav et al J 12009 ; IBouwens et~aHl2010bl ; 
iFinkelstein et alj|2010h 

The redshift evolution of the luminosity function of 
LAEs is different. The Lya luminosity functions of LAEs 
are observe d to be remarkably constant betwe en z = 3 and 
z = 6 (e.g. iHu et al]|l998l ; lOuchi et alj|2008h. after which 
the number density decreases at z > 6 (IKashikawa et al.l 
120061 : lOta et alj|20ba l20ld ; IKashikawa et al.ll201ll . but also 
see Hu et al. 2010). The first of these observations may be 
another consequence of the decreasing d ust content of star 
forming galaxies towards higher redshift (|Haves et alj|20l1l ; 
iBlanc et al.ll201lh . The second observation has received sig- 
nificantly more attention because it may signpost the exis- 
tence of large regions of intervening neutr al intergalactic gas, 
which are opaque to the Lya photons dHaiman fc Spaand 
Il999l ; iMalhotra fc Rhoadsl I2OO6I ; IKashikawa et all 120061 ). 



However, a proper interpretation of the redshift evolution 
of LAE luminosity functions at all redshifts requires un- 
derstanding the detailed radiative transfer of Lya pho- 
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Diikstra et al.ll201llh and ionized intergalactic medium 


(e.g. iDiikstra et al. 


2007b; McQuinn et al.ll2007; Ilicv et al. 


2008; 


Zheng et al. 


2010a; iDaval et al. 201ll; lLaursen et al. 


2011; 


Diikstra et al. l201ll, also see e.g. Fernandez & Ko- 



matsu 2008). 

LAEs are selected differently in theoretical models and 
observations. Observationally, LAEs are defined by their lo- 
cation in several 2-dimensional color-colour spaces. For ex- 
ample, Ouchi et al. (2008) define z — 3.1 LAEs by requiring 
that V - NB503 > 1.2 and [(V < V 2ct and B - V > 0.5) 
or (V U 2ct and B — Via > 0.5), where V 2a = 27.7 de- 
notes the 2a limiting magnitude of the V-band images] . The 
first requirement corresponds to having an excess flux in the 
narrowband filter denoted by NB503, while the broad-band 
selection corresponds to the requirement of having break 
in the spectrum blueward of the Lya resonance. The re- 
quirement that the object have a strong narrowband excess 
translates approximately to a minimum restframe equivalent 
width (REWmin) of the Lya line. In detail REW m in depends 
on the S /N at which an object is detected (see e.g. Gronwall 



et al. 2007 for an extended discussion on this and other com- 
plications regarding the flux dependence of the survey vol- 
ume). Representative values of REW m i n are often reported 
in papers discussing measurements of the LAE luminos- 
ity function, but these numbers should be interpreted with 
caution (e.g. Ouchi et al. 2008). However, when describing 
the LAE population in numerical or semi-analytic calcula- 
tions, theoretical models generally literally adopt the quoted 
minimum val ues of REW m j n to sel e ct LAEs from model 
galaxies (e .g. Le Delliou et al. _ 20061; Dijkstra et al.l 2007a : 



| 2010l; iNagamine et al 
l2010al: IShimizu et al 



Mao et all ^0071; iMcQuinn et'alj I2OO7J; iKqbavashi et al 



■20ld : IDaval et alj|201ll; IZheng et~. 
201 ll ; iForero- Romero et all l201ll ).~Ln 



addition, some theoretical papers only apply flux thresholds, 
which formally separates them further from observations. 

In this paper, we investigate the impact of simplified se- 
lection criteria (i.e. joint REW and luminosity cuts) on the 
'predicted' number density of LAEs. In particular, we gener- 
ate phenomenological models for the number density of star 
forming galaxies as a function of UV and Lya luminosity. We 
then investigate whether we can reproduce the LAE data if 
we define LAEs by a (redshift dependent) Lya luminosity 
and equivalent width (EW) threshold. Our empirical mod- 
els are taken from observed drop-out luminosity functions at 
3 ^ z ^ 6, combined with observational constraints on the 
prominence of Lya emission lines in drop-out galaxies. The 
Lya equivalent wi dth distributio n for L BGs has been mea- 
sured at z = 3 by IShaplev et alj ([2003), and at 3 < z < 7 
for a slightly smaller sample by Stark et al. (2010, 2011). 
We then investigate whether it is possible to reproduce the 
observed Lya luminosity functions (i.e. the number density 
of LAEs as a function of Lya luminosity and UV magni- 
tude) and equivalent width distributions, based solely on 
our knowledge of this LBG population. If this is not pos- 
sible, then the analysis implies that utilising simplified se- 
lection criteria can have a strong impact on the predicted 
number of LAEs. Indeed, our results indicate that as models 
of LAEs mature, selection criteria that resemble those ap- 
plied on the actual data will increasingly have to be taken 
into consideration. 

The outline of our paper is as follows: in §f2]we describe 
our formalism for 'predicting' the Lya luminosity function, 
and its uncertainties. In §[3] we present our results for Lya 
luminosity functions at z <J7, and provide a discussion of 
uncertainties in § [4] We also make predictions for LAE lu- 
minosity functions out to z = 10 in § [5] Finally, in § [6] we 
present our conclusions. The cosmological parameter val- 
ues used throughout o ur discussion are (Sim , SIa, Sib, h) = 
(0.27,0.73,0.046,0.70) jKomatsu et alj|2009h . 



2 EMPIRICAL RELATIONSHIP BETWEEN 
LBGS AND LAES 

The range of UV-magnitudes probed by samples of LBG and 
LAE^I currently overlap. Furthermore, recent Lya radiative 
transfer modeling shows that Lya photons esc ape anisotrop- 
ically from young, dusty, simulated galaxies (|Laursen et alj 
2009). This implies that the observed narrow-band excess 
of a star forming galaxy - and hence, whether it would be 



For many LAEs the UV continuum is also detected. 
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selected into a sample of LAEs - depends on its orienta- 
tion relative to the observer. Furthermore, modelling of the 
observed spectral profiles of LBGs near the Lya line sug- 
gests that LBGs observed not to have any Lya emission 
are not a separate class of galaxy. Rather, these are also 
intrinsically strong Lya emitters, in which radiative trans- 
fer effects tr ansform the Lya emission line into an absorp- 
tion feature ilSchaerer fc yerhammell2008l ; lAtek et al J 120091 ; 
iDessauges-Zavadskv et all 20100 Together these observa- 
tions indicate that the LBG and LAE populations are intrin- 
sically the same, and are only separated as a consequence of 
Lya radiative transfer through opaque (dusty) media. 

The possibility that Lya radiative transfer plays an im- 
portant role in separating LAEs from LBGs is underlined by 
the observation of low surface brightn e ss, extended Lya ha - 
los around LBGs l|Steidel et all 1201 ll ). Isteidel et ail i|201ll ) 
argue that the Lya radiation in these halos was emitted as 
nebular emission, which is then scattered to the observer in 
the (outflowing) circumgalactic medium. When the flux in 
these halos is properly accounted for, the total observed Lya 
flux places most LBGs in the LAE category. 

Thus, current evidence suggests that the LBG and LAE 
populations should be considered within the same theoretical 
framework. However, before this can be reliably pursued it is 
important to understand the observational selection criteria 
that define LAEs. 

2.1 The Formalism 

Our goal is to connect the observed numbers of LAEs and 
LBGs. In this section we therefore begin by introducing the 
formalism used in this paper. The number density of LAEs 
with Lya luminosities in the range L a ± dL a /2 is 

<&(L a )dL a = dL a x F x (1 + 5v) x 

/■Mmax 

/ dM uv <t>{M m )P{L a \M uv ), (1) 

where <f)(M uv )dM\jy denotes the number density of galax- 
ies with absolute AB-magnitude range M uv ± dMuv, and 
P(L a \M uv ) is the conditional probability density function 
(PDF) for L a for a given Muv- The parameter F is a nor- 
malisation factor which is discussed in § 12.31 The factor 
(1 + Sv) = Ay accounts for cosmic variance, and we as- 
sume that Ay is drawn from a Gaussian distribution with a 
standard deviation ay = 0.36 (estimated following the pro- 
cedure of Somerville et al. 2004, see Moster et al. 2011 for an 
update of this workjf). The integral is taken over the range 
M min = -30.0 to M max = -12.0. 

2 How frequently such a 'transformatio n' occurs can depend o n 
the metallicity of the interstellar gas (e.g. iFinkelstein et al.ll201ll) . 

3 The cosmic variance recipe given by Somerville et al (2004) 
formally applies to spherical volumes. The narrowband survey 
of Ouchi et al. (2008) probes volumes that are close to cubical, 
and therefore Somerville et al. (2004) should provide reasonable 
estimates for their cosmic variance. For the appropriate survey 
volume of 0.5 - 0.9 X 10 6 Mpc -3 , Figure 3 of Somerville et al. 
(2004) gives odm = 0.03 — 0.06. To get a conservatively large 
estimate for <r, we adopt a linear bias paramet er of b = 6 which 
is on the high end of the observed range (e.g. IShim asaku et alj 
120031 : iGawiser et~al"1l2007l; iKovac et aI1l2007l ; iGuaita et al.ll2010l) . 
This gives us a ~ 0.18 - 0.36. 
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Table 1. Adopted parameters for <p(M uv ). 



redshift 
z = dfl_ 
z = 4 b 
z=~ 


4>* (KT 3 cMpc" 3 ) 
1.7 ±0.5 

1.3 ±0.2 
1 4+0.6 


^fjv 
-21.0 ±0.1 

-21.0 ±0.1 
-20.2 ±0.2 


a 

-1.73 ±0.13 
-1.73 ±0.05 
-1.74 ±0.16 


2 = 7^ 

2 = 8£ 

2 = 9 C 


1 - i -0.7 
1 1 + 1 ' 7 

1 1 + 1 ' 7 
1 - 1 -0.7 


-19.8 ±0.4 
-19.45 
k, ~ 19 - 6 


-1.74 
-1.74 
-1.74 



a From Reddy & Steidel (2009). 
b From Bouwens et al. (2007). 

c From Bouwens et al. (2008). The slope a was kept fixed at 
the value that was inferred from the lower redshift observa- 
tions. To obtain the lower limit on M uv , the value for <j>* at 
2 = 9 was assumed to be the same as at redshift 7. 

From Bouwens et al. (2010a). The constraint on -M uv was 
obtained by assuming no evolution in <j>* and a. 

Table 2. Parameters related to detection thresholds in the 
narrowband surveys. 



redshift 


REW mil £ 


r , /erg 
a, mm/max \ s 


t 


2 = 3.1 


64 A 


2 / 50xl0 42 




2 = 3.7 


44 A 


4 / 40xl0 42 




2 = 5.7 


27 A 


2.5 / 40xl0 42 





a Taken from Table 3 of Ouchi et al. (2008). 
b The Lya luminosity functions of Ouchi et al. (2008) has 
bins with a width of 0.2 dex (in logioL a ). To estimate L aim i n 
(£max), we subtracted (added) 0.1 dex from (to) the Lya lu- 
minosity of the faintest (brightest) bin. 

2.1.1 The Drop-Out Luminosity Function 4>(M uv )dM\j\r 

We assume that the LBG luminosity function <^(Muv)dM uv 
is described by a Schechter function, the parameters of which 
(0*, Mu V , a) we take from the literature. Table [T] summa- 
rizes the redshift dependence of our adopted parameters, 
and the references from which these were taken. Note that 
formally, the Lya luminosity functions have been determined 
at 2 = 3.1, 2 = 3.7 and 2 = 5.7. We have interpolated the 
UV luminosity functions to these same redshifts by assum- 
ing that M T * TV evolve s as M uv = -21.02 + 0.36(2 - 3.8) 
(|Bouwens et al.l 120081 ') while keeping the other parameters 
fixed. 

2.1.2 The Conditional Probability P(L a \M uv ) 

The absolute AB UV-magnitude relates to the UV lu- 
minosity density Z/u v,^ (in erg s -1 Hz -1 ) as Muv = 
— 2.51ogI/uv,i/ + 51.6 (|Ouchi et aklbOOot ). Because Lya lumi- 
nosity is simply the product of rest frame equivalent width 
(REW) and luminosity density (La in erg s^ 1 A" 1 ), we can 
express Lya luminosity as a function of REW and Z/uv.i/ as 
L a = C x REW x Luv.k- The constant C = ^(^f)"^" 2 , 
in which u a = 2.47 x 10 15 Hz, X a = 1216 A, = 
dlog Z/A/dlog A (i.e. L\ oc X 13 and L v oc u^ 13 ^ 2 ), and Auv = 
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1700 A denotes the restframe wavelengt h at which the UV 
conti nuum flux density was measured (|Diikstra fc Westral 
2010). Throughout this work we assume /3 = —1.7 (see § 14,211 , 
We can thus recast the conditional probability P(L a \M uv ) 
as a function of the REW-PDF as 



dRE\V c dRE\Y v 



P(L a \M u 



P(x\Mvv) m:xY 




dh 



X £ (x 

min j ^max 

otherwise, 



(2) 



where x =REW=L a /[C x Luv,^], and x m i n denotes the 
minimum equivalent width REW m i n that a given narrow- 
band survey is sensitive to (for smaller values of REW 
the narrow band excess would be too small for the galaxy 
to make it into the sample of LAEs). Numerical values 
of 

^min =REW m in that approximately represent the color- 
cuts adopted by Ouchi et al. (2008, see their Table 23) are 
given in Table [2] For the maximum REW we have assumed 
a; max =REWmax = 300 A, but note that our results do not 
depend precisely on this number (see § I4.2|l . 

The conditional PDF for the REW [P(REW|M UV )] has 
been measured for LBGs. In Appendix lAll we show that the 
Lya REW PDF can be well described by an exponential 
whose scale length depends on Muv and z 



P(REW|M UV ) = JVexp | 



-REW 



(3) 



VREW C (M UV ) 

Here j\f denotes a normalization constant, which we choose 
so that all drop-out galaxies have — ai ^REW^REW max 
(see Appendix Al). In this expression the factor ai = 20 
A for Muv < -21.5, and oi = 20 - 6(A/ UV + 21. 5) 2 A for 
—21.5 Muv ^ —19.0, and we freeze the evolution of a\ 
for Muv > -19.0. Furthermore, REW C (M UV ) = REW c ,o + 
^AM + ^§^Az with AM = M vv + 21.9, and Az = 
z - 4.0. The best fit values are REW c ,o = 22 ± 3 A, and 
= 6 ± 4 A. Equation (|3j) ensures that the fraction 



dREW c 
dM uv 



of drop-out galaxies with Lya in emission (i.e. REW .> 0) 
depends on Muv, as is observed. We assume throughout that 
this fitting formula applies only in the observed range of UV 
magnitudes, and also freeze the evolution of REW c (Muv) for 
Muv > —19.0 (see § I4.2l for a discussion of the uncertainties 
this may introduce) . The choice of the functional form that 
approximates the data is quite arbitrary. In Appendix IA2 
we investigate an alternative parametrization, and find that 
our results are not significantly affected. 



2.2 Comparison to Data 

We next compute 'predicted' Lya luminosity functions by 
combining Eq Q] Eq H and Eq [3] at z = 3.1, z = 3.7 
and z = 5.7. To facilitate the comparison with the data, 
we compute the quantity ^{L a )d\ogL a which denotes the 
number density of LAEs in the range logL a ± (dlogL a )/2. 
The units of ^(L a ) are cMpc -3 [log LJ -1 , and we have 
*(L a ) = \nWL a ®{L a ). 

Each model is described by three parameters that quan- 
tify the Lya REW distribution (REW c ,o, d j^ c , dR ^ c ), 
the normalization parameter F, and three Schechter 
parameters plus one cosmic variance parameter at 
each redshift bin. We therefore need to explore a 16- 
dimensional parameter space (3 + 1 + 3x4 = 16), 
which we do with a Monte Carlo Markov Chain 
(MCMC) method. We characterize each model by 



the parameter vector P = (REW c ,o, —— — • — — 

, F,ct3,..,6,4>t 6,-Muv,3,...,6 ! <5v,3,...,6)- Here, we have 

adopted a notation where «3 denotes the value of a at 
z — 3. We compute the posterior probability for each model 
as P(P) oc £[P]P(P), where £[P] = exp[-0.5x 2 ] denotes 
the likelihood, in which \ 2 = ^^ Vdata (modeh — datai) 2 /u 2 . 
The function P(P) denote s the prior PDF for these 
parameters (e.g. ICowanlll997h . 

We have assumed that P(P) is a multivariate Gaussian, 

i.e. P(P) = jVexp [ - |(P - ^p) T C- 1 (P - M p)] , where N 
denotes the normalization factor. The vector /ip contains 
the best fit values for each of the parameters (e.g. from Ta- 
ble 1 we have /i a ,3 = —1.73, and in § 12.1.21 we find that 
£tREW,c,o = 22 A). Note that we do not assume any prior 
knowledge of F or . The covariance matrix C con- 

tains the measured uncertainties on the parameter^]. Table Q] 
summarizes the assumed redshift evolution for M{j v , cj>* and 
a and their uncertainties. We then compute marginalized 
PDFs for the parameters F and dr ^ w by marginalizing over 
the other 15 parameters. The data that we use for the fits is 
from Ouchi et al. (2008). 

Finally, we summarize other model parameters that we 
do not vary as part of the MCMC calculations in Table [5] 
These parameters include: (i) the approximate minimum 
equivalent width (REW m i n ), and (ii) the minimum Lya lu- 
minosity (L a ) to which the narrowband survey of Ouchi et 
al. (2008) was sensitive. 



2.3 The 'normalization' parameter F 

The 'normalization' parameter F scales the predicted lumi- 
nosity function up and down, and so can be interpreted as 
the ratio of the observed to predicted number density of 
LAEs. In this paper we focus on the value of F that ar- 
rises from fitting to LAE and LBG luminosity functions, as 
a means of testing the validity of using a simple cut in Lya 
luminosity and REW to represent the real observational se- 
lection criteria. In a case where the simple cuts accurately 
represent the true selection criteria (and hence models pro- 
vide a faithful representation of the data) we therefore ex- 
pect that F — 1. However any deviation of F from unity 
indicates that the simple cuts are not providing an adequate 
description of the observational selection criteria. 



4 The covariance matrix in this case is a 16 X 16 matrix whose en- 
tries are given by dj = (JiOjPij. Here Ui denotes the uncertainty 
on parameter '«', and pij denotes the correlation coefficient be- 
tween parameter i and j. While these correlation coefficients are 
generally not given, the constraints on the parameters a, Mfjy , <f>* 
at a given redshift are strongly correlated. We assumed through- 
out that p a ,M* = Pm* ,£* = Pa,<t>* = 0.9 at each redshift. We 
found that this decently reproduces the shape of the 68% and 
95% likelihood contours for different Schechter parameter combi- 
nations as given by Bouwens et al. (2007). We assume that the 
other parameters are not correlated, i.e. pij = 0. By definition 
pa = 1 for all V. 
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a 




Lycx luminosity (erg/s) 




Figure 1. We plot Lyc* luminosity functions at z = 3.1 (upper left), z = 3.7 (middle) and 2 = 5.7 (upper right). The soiid ftZacA; Zines 
show our best fit model, which we calculated following the procedure described in § 12.1.11 -§ 12.21 The data of Ouchi et al. (2008) are 
indicated as blue circles in all panels. For completeness, we have also shown some other data sets, which were not included in our fits 
(see text). The lower left panel show marginalized posterior PDF for the parameter F. The lower right panel shows the 68% and 95% 
contours in the F -5 — - plane, which we obtained by marginalizing over all other parameters. 
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Figure 2. The predicted UV luminosity function of LAEs for our best-fit model (black solid lines), compared to the data of Ouchi et al. 
(2008). The agreement between our model and the data is excellent at z = 3.1 (left panel, note that this is not a fit), and reasonable at 
z = 3.7 (central panels). Our model underprcdicts the UV luminosity function at z = 5.7, which may be related a slight overabundance of 
large REW (REW J>200 A) systems in our model (see text). For completeness, we have also shown the best-fit drop-out galaxy luminosity 
function (black dashed lines), and the best fit Schechter function of the UV luminosity function of LAEs as derived by Ouchi et al. (2008, 
red dotted lines). 



3 RESULTS 

3.1 The Lya Luminosity Functions 

In Figure [T] we show Lya luminosity functions at z = 3.1 
(upper left), z = 3.7 (middle) and z = 5.7 (upper right). The 
solid black lines show our best fit model (see below), which 
we calculated following the procedure described in § 12.1.11 



-§ 12.21 The data of Ouchi et al. (2008) is indicated as blue 
circles in all panels. For completeness, we have also shown 
the z = 3 data from Rauch et al. (2008, green circles) and 
Gronwall et al. (2007, red circles), and the z = 3.7 data from 
Dawson et al. (2008, red squares). These other data points 
were not included in our fits. 

The best fit model is in excellent agreement with the 
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Figure 3. The model equivalent width PDF (black thick solid lines), compared to the best estimate of the observed REW PDF at 
z = 3.1, z = 3.7 and z = 5.7 (Ouchi et al. 2008), for all their photometrically selected LAEs (red filled squares), and for spectroscopically 
identified LAEs (blue filled circles). The vertical dotted lines indicate the minimum REW of the Lya emission line that galaxies need to 
have in order to be selected as LAEs by Ouchi et al. (2008). 



data, and is described by the parameters REW C o = 23 
A, = 7 A ' = 6 A, F = 0.53, a 3 = -1.65, 

(j>3 = 1.9 x 1(T 3 cMpc- 3 ,Muv,3 = -20.9, <5 v ,3=+0.04; a 4 =- 
1.70, cj>l = 1.3 x 10" 3 cMpc -3 ', M; V]4 = -20.9, <S V ,4=-0.10, 
a 6 = -1.73, 4>l = 1.5 x 10 -3 cMpc" 3 , M uv , 6 = -20.4, 
Sv,6 = 0.01. It is interesting - though likely a coincidence 
- that this model is in reasonable agreement with the data 
of Rauch et al. (2008, green circles). These authors detected 
27 ultrafaint Lya emitters in a long-slit, 92 hrs, observa- 
tion with the ESO VLT FORS2 spectrograph, and probed 
an effective volume of ~ 10 3 cMpc 3 . In such a small volume 
cosmic variance i s sign ificant. Following the prescription of 
ISomerville et al] |2004) for estimating cosmic variance, we 
find that the uncertainty on the number of detected galax- 
ies is on ~ 0.8(N). These uncertainties are denoted as grey 
dashed lines because this estimate is likely not accurate for 
the highly elongated survey volume probed by Rauch et al. 
(2008, see Munoz et al. 2010 for a discussion of cosmic vari- 
ance in pencil beam surveys). 



The lower left panel shows the marginalized posterior 
PDF for the parameter F. Based on available data for the 
drop-out galaxy population we find that we overpredict the 
number density of LAEs, and so need to multiply our pre- 
dictions by a constant factor of F ~ 0.43 to match the data. 
Formally, F = 1 is ruled out at ~ 95% CL. As noted previ- 
ously, in a case where sharp REW and luminosity cuts among 
the population of drop-out galaxies can be used to reproduce 
the observed luminosity functions of LAEs we expect F = 1. 
The fact that the data are inconsistent with F = 1 therefore 
demonstrates that the use of sharp luminosity and REW 
cutoffs does not adequately describe the observational selec- 
tion of LAEs. This represents the primary conclusion of our 
paper. 



In the lower right panel we plot the 68% and 95% 
contours in the F — dR ^ ; J Vc plane, which we obtained by 
marginalizing over all other parameters. The data clearly 
prefer models with dR ^J Vc > 0, in excellent agreement with 
the findings of Stark etal. (2010, 2011, also see Fig©. 



3.2 The UV Luminosity Functions 

We next take our best-fit model and compare it to other 
observed properties of LAEs. First, we 'predict' the UV lu- 
minosity function of LAEs, < I > lae(wiuv, z) 

$lae(™uv, z) = <j>(M uv , z) x J(M UV , z), (4) 

where <f)(M uv ,z) denotes the LBG UV luminosity function 
(introduced in §[5J, and J(Muv, z) denotes the fraction of 
drop-out galaxies that have REW^ REW m in and a total 
Lya flux greater than L a ,min (see Table [2] for numerical val- 
ues of L a , min and REW m in). Note that this latter constraint 
is important. For fainter drop-out galaxies, the requirement 
that the Ly-a flux be sufficiently large can translate to a 
higher required minimum REW (also see Zheng et al. 2010a). 
The fraction J-"(Afuv, z) is then obtained by integrating over 
the REW-PDF. 

Figure [5] shows the best-fit drop-out galaxy luminosity 
function (black dashed lines), the best fit Schechter function 
of the UV luminosity function of LAEs as derived by Ouchi 
et al. (2008, red dotted lines), and our predicted UV luminos- 
ity function of LAEs for the best-fit model shown in Figure [T] 
(black solid lines). The 'breaks' in the UV-luminosity func- 
tions at Muv > — 19 arise because Lya REWs greater than 
our quoted REWmin are required to render the Lya flux 
large enough (i.e. L a ^ L a , m in, see Table [2]) to be detected 
from fainter UV magnitudes. For example, in order for an 
object to be detected in Lya from an Muv = —18.5 galaxy, 
we need REW> 100 A, and these objects are significantly 
rarer. Note that the 'sharpness' of the break increases to- 
wards lower redshift. This is because the scale-length of the 
REW-PDF increases with redshift (see Fig 1), and objects 
with REW> 100 A are rarer at lower redshift. 

The agreement between our model and the data is ex- 
cellent at z = 3.1 (left panel, note that this is not a fit), 
and reasonable at z = 3.7 (central panels). In the z = 3.7 
case, the observed UV-LF suffers from incompleteness at 
Muv > —20.0, and so our formalism therefore overpredicts 
only the brightest UV point. Our best-fit z = 5.7 model 
clearly underpredicts the LAE UV-LF at z = 5.7 (right 
panel) at —20 < Muv < —21.4. The reason for this is not 
clear. It is possibly related to the fact that we overproduce 
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the number of large REW (REW ^200 A) systems (see be- 
low), although this discrepancy in the REW distribution is 
quite small. 

3.3 The Lya Equivalent Width Distribution 

We also use our formalism to predict the REW distribution 
(derived in Appendix [B} for our best-fit model 

P(REW, z) = 

M j a ' maX p(REW|M U v,c,2)0(M uv ,c^)dlog lo L Q , (5) 

where M denotes a normalization constant, <f)(AI uv , z) again 
denotes the LBG UV luminosity function (introduced in §|5), 
and P(REW|M uv , c ,z) denotes the Lya-REW PDF that is 
observed for LBGs. The relation between Lya luminosity, 
UV continuum flux density, and REW uniquely determines 
the absolute UV magnitude Muv.c at fixed L a and REW. 
This equation states that at a given Lya luminosity, the 
probability of observing a galaxy with Lya REW is the sum 
of all possible Muv, weighted by their number density. 

Figure O shows our model equivalent width PDF [black 
thick solid lines). The red filled squares show the observed 
REW PDF for all photometrically selected LAEs (also shown 
as the black histograms in Fig 23 of Ouchi et al. 2008). The 
blue dashed histograms show the observed REW distribution 
of spectroscopically confirmed LAEs. At z = 3.1 (left panel) 
our model underpredicts the observed number of drop-out 
galaxies with REW> 175 A. The agreement between our 
model and the data at z = 3.7 is excellent. As eluded to pre- 
viously, our model slightly overpredicts the number of large 
REW systems (REW> 200 A) at z = 5.7, which partially 
explains why we underpredict the observed UV-LF of LAEs 
at z = 5.7. 



4 DISCUSSION 

In the previous section we showed that the observed 
REW distribution of LAEs can be reproduced well out to 
REW=200 A using sharp cuts in REW and Lya luminos- 
ity among the drop-out population. However our empirical 
procedure leads to LAE luminosity functions (both UV and 
Lya) that are overpredicted by a factor of 1/F ~ 2.5. As 
already noted, the difference of F from unity indicates that 
simple selection cuts in luminosity and REW do not ade- 
quately represent the LAE selection. The goal of this dis- 
cussion is to explore the observational biases and model as- 
sumptions that may cause F to be less than 1. 

4.1 Discussion of Observational Biases 

Firstly, we note that using the drop-out galaxy population 
to constrain the number density of star forming galaxies 
as a function of their UV and Lya luminosity, does not 
miss the small fraction of 'red' star forming galaxie s (e.g. 
those with UV slopes ft J> — 0.5. lBouwens et afll2009l ). since 
the drop-out galaxy luminosity functions u sed in this pa- 
per h ave been corrected for this bias (see iBouwens et al.l 
I2009T ). Secondly, narrowband surveys can pick up galaxies 



such as ULIRGs which do n ot make it into drop-out sur- 
veys l|Nilsson fc M0lleri l201lT ). However, the ULIRG frac- 
tion among LAEs drop s dramatically to <10% at z > 2.5 
(jNilsson fc M0llerll2Ol j ). and this is unlikely to be a signif- 
icant effect. A hypothetical population of LAEs whose Lya 
emission is powered predominantly by gravitational heating 
(as in e.g. Birnboim & Dekel 2003, Dijkstra 2009, Dayal et 
al. 2010), would also contribute to LAE samples. However 
these known potential observational biases would enhance 
the true number density of LAEs, and thus lead to a value 
of F > 1. 

Ouchi et al. (2008) discuss possible explanations why 
F > 1 (see § 14.31 for a more detailed comparison with Ouchi 
et al. 2008). Some of these possibilities could also result in 
F < 1. These include: (i) Systematic uncertainties in drop- 
out galaxy luminosity functions, based the observed scatter 
in the z = 6 drop-out LFs obtained by different groups, (ii) 
The fact that drop-out LFs at a particular redshift are mea- 
sured over a much broader redshift interval (Az ~ 1) than is 
probed by narrowband surveys. These differences may intro- 
duce extra uncertainties when comparing LAE and drop-out 
galaxy populations. 

We have investigated whether our results are dominated 
by observations at a particular redshift, and repeated our 
analysis based on the LAE luminosity function in individual 
redshift bins (while fixing dRE JJ crit = 6 A). We find that 
F = 0.34+° :tr at z = 3.1, F = 0.44±g;£? at z = 3.7 and 
F = 0.52lg23 at z = 5.7, where the uncertainty on F is 
dominated by cosmic variance. Thus, our result is not domi- 
nated by any redshift bin. The systematic effect that causes 
F < 1 appears to operate at all redshifts. 



4.2 Discussion of Model Assumptions 

In this section we discuss the effects of different assumptions 
on our results, with particular focus on the conclusion that 
F< 1. 

Assumption 1. In our fiducial model we assume that 
the fitting formula for REW (Eq[3]) applies only in the ob- 
served range of UV magnitudes, and freeze its evolution 
for Afuv > -19.0. On the other hand, Stark et al. (2010) 
find that the observed evolution in the REW-PDF continues 
down to Muv = —18.5 (see the left panel of their Fig 13), 
albeit with large uncertainties. Had we extrapolated our 
fitting function down to fainter UV magnitudes, then we 
would have allowed more UV-faint, large REW galaxies into 
our sample. However, if we include these UV-faint sources, 
then we would push our constraints on F to lower values, 
which would rule out F = 1 at greater significance. Con- 
versely, the uncertainties on the REW-PDF are large at 
Muv 55= —19.0. If we had frozen the evolution of the REW- 
PDF at Muv ^ —19.25, then we would have found a larger 
value of F, especially when combined with a truncation of 
the REW-PDF at REW> 150 A (see below). 

Assumption 2. In our fiducial model we chose 
■l max =REW max — 300 A to be the maximum possible Lya 
REW. This value is close to the largest REW in the sam- 
ple of Ouchi et al. (2008). Theoretically, the Lya REW can 
reach ~ 1500 — 3000 A for metal free galaxies forming stars 
with a top-heavy IMF (Schaerer 2003, Johnson et al. 2009, 
Raiter et al. 2010), or cooling clouds (Dijkstra, 2009), and 
could be boosted to even larger values if dust preferentially 
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supp resses the UV continuum (|Neufeldlll99ll :l Hansen fc Oh] 
I2006T ). Adopting larger values for REW max would boost the 
overall Lya emissivity of star forming galaxies, which would 
again reduce our best-fit values for F. 

We note that there is limited data to support the as- 
sumed exponential form of the REW-PDF at REW ^150 A. 
Indeed, among the ~ 800 LBGs in the sample of Shapley et 
al. (2003), only 4 (1) have REWJs 150 (175) A, so that the 
observed REW-PDF among LBGs is very uncertain at these 
values. If we truncate the REW-PDF at REW max = 150 A, 
then we find F = 0.64±o*°, and F = 1 is onl y mled out 
at ~ 89% CL. If we further combine this with 'freezing' the 
evolution of the REW-PDF at M UV > -19.25, then we find 
F = 0.70±g;f|, and F = 1 is only ruled out at - 77% CL. 
While these modified assumptions do not fully resolve the is- 
sue, they do illustrate that our finding of F <C 1 depends on 
the uncertain REW-PDF at REW ^150 A and faint M uv . 
We note that if we assume REW max = 150 A, then our 
model does not produce LAEs with REW^ 150^4, although 
these objects are observed. However, considering just LAEs 
that have been confirmed spectroscopically, this appears to 
be a significant problem only at 2 = 5.7, where the uncer- 
tainties on measured REWs are large (see Table 2 of Ouchi 
et al. 2008). 

Assu mption 3. In our fiducial model we assumed that 
P = -1.7. IStark et all |2010l ) found that /3 = -1.6 for drop 
out galaxies with REW< 50 A, and f) = -2.0 when REW^ 
50 A for -21.5 < Muv < —20.5. For fainter galaxies with 
REW< 50 A, /3 approaches —2.0. Our cho i ce fo r /3 ma y 
be slightly too high (also see iBouwens et al.1 [20091 . l2010bl ). 
However, if we decrease j3 then we increase the total Lya 
flux for a fixed REW and Muv, and thus the overall Lya 
emissivity of the drop-out galaxy population. Decreasing f3 
would therefore again lower our overall best-fit value for F. 

4.3 Comparison to Previous Work 

Our study bares similarities to that of Malhotra & Rhoads 
(2002). These authors also compared observations with 'pre- 
dicted' LAE number counts and REW distributions, which 
they obtained by combining existing constraints on the 
faint-end of the 2 ~ 4 drop-out galaxy luminosity function 
with a theoretical model for the Lya REW-distribution. 
However, modeling the observed equivalent width of the 
Lya emission line is a complicated task which depends on 
detailed radiative transfer of Lya photons through both 
ISM and IGM. Our calculations completely circumvent 
this complication by utilizing empirical distributions, which 
represents an important i mprovement. In agre e ment with 
our findings, the study of iMalhotra fc Rhoadsl (I2002T ) also 
overpredicted the number density of LAEs. Their offset 
was by an even larger factor of ~ 6 — 12, implying that 
only ~ 7 — 15% of all LBGs need to be LAEs. However 
our formalism already includes the observation that only a 
fraction of the drop-out population has strong enough Lya 
emission to qualify as a LAE. Thus the concept of separate 
galaxy populations cannot be invoked to explain the offset 

Ouchi et al. (2008) found that the UV-LF of the LAEs 
at 2 = 5.7 lies remarkably close to the UV-LF of drop- 
out galaxies at 2 = 6, which suggests that ~ 50 — 100% 



of the drop-out galaxies at z=6 would qualify as LAEs 
(also see Shimasaku et al. 2006). Phrased alternatively, the 
UV-LF of LAEs implies that - 50 - 100% of the drop- 
out galaxies have a Lya emission line whose REW ex- 
ceeds REWcrit ~ 20 A. This conflicts with spectroscopic 
observations of LBGs at that redshif t which suggest the 
observed number is closer to ~ 30% dShaplev et al.l [20031 : 
Sta nwav et al. I l2007l ; iDow-Hygelund et al.ll2007r ). 

The 2 = 5.7 UV-LF of Ouchi et al. (2008) therefore 
implies that F ~ [50 — 100%]/30% = 1.6 - 3.3, i.e. they find 
more LAEs than expected from drop-out galaxy populations. 
At first glance this conflicts with our finding that F < 1. 
However, our constraint on F was derived by considering 
the Lya LFs at 2 = 3.1, 2 = 3.7 and 2 = 5.7. We showed 
in § 14.11 that F = 1 was only excluded at the 1 — a level if 
we had only considered the 2 = 5.7 Lya LF. Importantly, 
the z = 5.7 UV-LF is the only luminosity function of the six 
LFs that we modelled for which F = 1 would have given a 
good fit (Figure [2] shows that our best-fit model undershoots 
the UV-LF by a factor of ~ 2, despite the fact that our 
model reproduces the Lya LF and REW distributions at 
this redshift as well as it does at other redshifts). 

The z = 5.7 UV-LF alone therefore appears to be con- 
sistent with direct spectroscopic observations of LBGs. This 
still seems at odds with Ouchi et al. (2008) who were con- 
cerned with explaining why F significantly exceeded unity. 
The discussion in Ouchi et al. (2008) used observational con- 
straints on the REW-PDF of the drop-out galaxy population 
that were available at that time. However recent data which 
contains larger samples of galaxies shows that the fraction of 
drop-out galaxies with REW> 20 A increases dramatically 
at Muv >-20.0, which is relevant when comparing to LAEs. 
If one accounts for this increase, then the discrepancy noted 
by Ouchi et al. (2008) becomes less serious. Indeed, in the 
recent compilation by Ono et al. (2011), the fraction of faint 
drop-out galaxies for which REW^ 25 A is ~ 55 ± 15%. It 
therefore seems likely that both this work and that of Ouchi 
et al. (2008) would conclude that the UV-LF of 2 = 5.7 
LAEs is consistent with direct spectroscopic observations of 
LBG. In contrast with previous studies, our constraints are 
derived from a model that uses more available data on the 
Lya REW-PDF observed in drop-out galaxies, considers the 
Lya luminosity functions, and studies different redshifts. 



5 EXTRAPOLATING TO REDSHIFTS Z ~ 7 - 9 

There are existing constraints on the LBG luminosity func- 
tions at 2 > 6. It is therefore interesting to take our best-fit 
model and predict the expected number density of LAEs at 
redshifts beyond those currently observed. Before doing so, 
we first point out that extrapolating the REW-PDF of our 
best-fit model to 2 > 6 leads to inconsistencies with some 
available data. 

The observed fraction of drop-out galaxies having Lya 
emission lines of REW^ x has been measured as a function 
of redshift out to 2 ~ 7 for galaxies with —20.5 < Muv < 
-19.5 for x = 75 A l|Stark et alj|2010h . and for galaxies with 
-20.25 < M uv < -18-75 for x = 25 A (|Pentericci et al.l 
l201ll ; ISchenker et al.|[2oTll ; IOno et al.ll201ll . also see Vanzella 
et al. 2011). In Figure [4] we compare the data from Stark 
et al. (2010, blue filled squares), and the data compiled by 
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Figure 5. Left Panel: Predicted number density of LAEs (in cMpc - 3 log _1 L a ) at z = 7 (solid red line). Right panel: Same as the left 
panel, but at z = 9. Within the range of luminosities that will likely be probed by future surveys, L a ~ 10 41 — 10 43 erg s — 1 , the agreement 
between previous model predictions (here dashed line is the GALFORM model from Nilsson et al. 2007, and the dotted line is the model 
of Dijkstra ct al. 2007a) and those inferred from the z=9 LBG population is quite good (see text). 



o 

«0.5 



• l REW>25 A; '-20.35 , <M U y<-lB.75;' data 
■ REW>75 A; -20.5<M UV <-19.5; data 
best fit model (-19.75) 




Figure 4. The fraction of drop-out galaxies with Lya emission 
lines with REW> function of redshift. Blue filled squares 

show the data from Stark et al. (2010) for galaxies with —20.5 < 
M;jy < —19.5 and x = 7 5 A, while red fill e d circ le s show the data 
from Stark et al] teOld) ; iPentericci et al] ll201ll1 ; ISchenker et al] 
j20111) ; |Pentericci et ahl (l201lf> compiled by Ono et al. (2011) for 
-20.25 < M uv < -18.75 and x = 25 A. The solid lines shows 
our best fit model (this model was fitted to the Lya luminosity 
functions, not this data), which provides and excellent fit to the 
data for our purposes (see text). The dotted lines show our best- fit 
model extrapolated to z ^ 6. The observed fraction at z = 7 falls 
well below this extrapolation. 



Ono et al. (2011, red filled circles) with our best-fit model 
of §[3] (solid lines). Our best fit-model slightly overpredicts 
(by a factor of ~ 1.2) the fraction of drop-out galaxies for 
which REW^ 25 A at z ^ 4. However, at these redshifts 
REW m in ^ 44 A, and our results depend weakly on this 
discrepancy. For our purposes, the agreement between our 
model and the data is excellent out to z ~ 6. We stress that 
this best-fit model was fitted to the Lya luminosity func- 
tions, and not these particular data points. This indicates 
that our model reproduces the observed redshift evolution of 
LAEs well, despite the fact that the overall predicted num- 
ber density is off by a factor of F . However observations 
indicate that a sudden drop occurs in the 'LAE fraction' at 



z > 6, and the linear extrapolation of our model does not 
reproduce this evolution. 

To remain consistent with the observed drop in the 

"' "'^Az 



'LAE fraction', we replace the term REW c ,o + — 
(yielding REW Cj 7 = 43 A ), in our best fit model with 
REWc.t = 8 A . This ensures that we reproduce the observed 
drop in the LAE fraction. Otherwise, we take the Schechter 
function parameters given in Table [TJ and the model pa- 
rameters from our best-fit model. In the left panel of Fig- 
ure [5] the red solid line shows the predicted number density 
of LAEs (in cMpc -3 \og~ 1 L a ) at z = 7, where we further 
assumed REWmin = 27 A(to facilitate the comparison with 
the z — 5.7 data by Ouchi et al. 2008). Our model pre- 
dicts the cumulative number density of LAEs brighter than 
L a = 10 43 ergs" 1 to be n(L a > 10 43 erg/s) ~ 10" 6 cMpc" 3 , 
which is below the observational constraints by Ota et al. 
(2010). However the observational constraints are uncertain 
given that their observed cumulative luminosity function is 
derived from only three objects. Our predicted luminosity 
functions are in turn affected by uncertainties in our model 
parameter vector P. For example, taking REW C ,7 = 15 A— 
which gives a Lya fraction among z ~ 7 drop-out galaxies 
of 35%- results in n(L a > 10 43 erg/s) ~5x 10" 6 cMpc" 3 , 
which is within 1 — ct of the best-fit value derived by Ota 
et al. (2010). Note that similar uncertainties will apply to 
any theoretical model, as these are likely to be -just like our 
parameter vector P-calibrated by lower-redshift data. 

The right panel of Figure [5] shows our predictions 
for redshift z — 9. We assumed Mf} v = —19-1, which 
we obtained by extrapolating the observed redshift evo- 
lution at low er redshift (z < 6) Mf} v = -21.02 + 
0.36(2: — 3.8) (|Bouwens et al]|2008l ). and we assumed that 
REWc,9 =REW C ,7. This extrapolated value is consistent 
with the derived low er limit Mf} v J> — 19.6 derived by 
iBouwens et al] (|2008l ) based on the absence in J-dropouts 
in the Hubble Ultra Deep Field (see Oesch et al. 2011 for 
constraints on the z = 10 drop-out galaxy luminosity func- 
tion using deeper and wider WFC3-IR data). For complete- 
ness, we have overplotted some theoretical predictions for 
the number density of LAEs at z = 9. The dotted line 
shows the predicted number density of LAEs at z = 8.8 by 
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Nilsson ct al. I J2007h who used the semi-analytic model GAL- 
FORM (|Cole et al.ll200oh . The dashed line was obtained from 
a simpler model in which some fraction /„ o f all baryons is 
conver ted into stars over a timescale CDcthub (|Diikstra et alj 
l2007aT ). These simpler models easily reproduce the observed 
number density of LAEs at z ^ 5.7, and have been used 
frequently in the recent literature. The theoretical predic- 
tions also agree well with each other, which is probably be- 
cause the models are calibrated by the same lower redshift 
data. The theoretically predicted luminosity functions de- 
crease more steeply than those obtained from the LBG pop- 
ulation. However, given the present-day uncertainties on the 
z = 9 LBG luminosity function, it is not clear how significant 
this difference is. For example, a steeper faint end slope of the 
z = 9 luminosity function would would reduce the discrep- 
ancy. Observations indicate that this faint end slope may 
indee d become steeper at higher redshifts (|Bouwens et al.l 
2011, also see Jaacks et al. 2011). 



6 CONCLUSIONS 

In this paper we have investigated the implications of the as- 
sumption often used in theoretical modelling that the LAEs 
are selected using cuts in REW and luminosity. These cuts 
are only coarse approximations to the detailed criteria that 
are employed in observational studies. To quantify the im- 
portance of the approximation, we investigated whether we 
can reproduce LAE data if we define LAEs using a (redshift- 
dependent) Lya luminosity and EW threshold, and a empir- 
ical model for the number density of star forming galaxies as 
a function of their UV and Lya luminosity. We constructed 
this model by combining observed luminosity functions of 
drop-out galaxies, with the observed rest-frame equivalent 
width (REW) probability distribution function (PDF) of 
drop-out galaxies at z = 3 — 7, and 'predicted' the resulting 
Lya luminosity function at z = 3.1, z = 3.7 and z = 5.7. We 
also use our formalism to predict the UV LF of Ly-a emit- 
ters, and the EW observed distribution of Ly-a emitters. 

As part of our analysis we demonstrate that the ob- 
served REW-PDF of z — 3 LBGs is well described by 
an exponential function at REW^ 0, i.e. P(REW) or 
exp[-REW/REW c (Afuv)]. The scale length depends on ab- 
solute UV-magnitude, and we use recent data from Stark et 
al. (2010, 2011) to constrain this dependence. Using this em- 
pirical distribution of REW in LBGs, we find that we can 
reproduce the observed REW distribution for LAEs quite 
well out to REW~ 200 A. However, in order to reproduce 
the LAE luminosity functions, we find that we must re-scale 
the predicted luminosity functions downward by a factor of 
F — 0.43to'.o7- Formally, a value of F = 1, which is ex- 
pected if the simple REW and luminosity cuts are accurate, 
is ruled out at ~ 95%. We found that this discrepancy can be 
reduced if we truncate the Lya REW-PDF at REW^ 150 
A, and 'freeze' its evolution at Muv J> — 19.25. The sam- 
ple of Shapley et al (2003) only contains 4 drop-out galaxies 
(out of 797) with REW^ 150 A, and observationally the 
REW-PDF is constrained very poorly at these large values 
for REW. For this truncated REW-PDF, which stops evolv- 
ing at M UV £ - 19.25, we find F = 0.70± ll, and F — 1 is 
ruled out only at ~ 77%. 

On the other hand, we found that the overall redshift 



evolution of the LAEs was reproduced very well by our em- 
pirical model. Encouraged by this result, we combine the 
best-fit model in our formalism with recent observed con- 
straints of the 'LAE fraction' and predict number densities 
of LAEs at z = 7 and z = 9. Current measurements of the 
LBG luminosity functions translate to Lya luminosity func- 
tions that are consistent with observed number counts of 
LAEs at z = 7, but which are still very uncertain. Neverthe- 
less the UV LF of LBGs can be used to provide empirical 
guidance for future surveys aiming to discover Ly-a galaxies 
at the highest redshifts. 

We conclude that modeling LAEs with simple REW and 
luminosity cuts can lead to (significant) changes to the pre- 
dicted number density of LAEs. Theorists will therefore need 
to make more careful account of observational selection in 
order to produce reliable models of the observed popula- 
tion. This will have to include application of the proper 
filter transmission curves to generate mock data from the 
models, and then correct for these filter transmission curves 
following the same procedures as followed by the observers 
whose data one tries to reproduce. We are attempting to 
address these issues in more detail in on-going work. Fur- 
thermore, our work has shown that to in order to address 
this issue in more detail, it will be important to reduce the 
observational uncertainties associated with the REW-PDF 
of drop-out galaxies at large REW and/or faint Muv- 
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Figure Al. The observed Lya rest frame equivalent width 
(REW) distribution of ~ 800 z ~ 3 LBGs of the sample of Shap- 
ley et al. (2003) is shown as the black histogram. This sample 
had a median Afyv = —22.3. This figure shows that the observed 
REW-PDF can be described well by an exponential function. 
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APPENDIX A: ACCURACY OF FITTING 
FORMULA FOR OBSERVED EW-PDF 

Al Fiducial REW-PDF Parameterization 

In our paper we approximate the observed Lya REW distri- 
bution using the functional form given by Eq(3] This func- 
tional form is motivated by several factors. 

• The observed REW-PDF of z = 3 LBGs is well de- 
scribed by an exponential function. This is illustrated in 
Fig IA11 where the histogram shows the observed number 
of LBGs as a function of Lya REW (Shapley et al. 2003). 
This sample had a median Mtjv = —22.3. The inset of this 
Figure shows that the best-fit scale-length associated with 
this exponential distribution is REW C ~ 29 A, for which the 
function is overplotted as the dotted line. 

• The scalelength REW C is observed to be a function of 
absolute UV magnitude, Mtjv (Shapley et al. 2003). This is 
illustrated in Figures[A2]and|A3l where show the REW-PDF 
for subsamples of the brightest ~ 400 LBGs (Fig |A2|) . and 
of the faintest 400 LBGs (Fig IX3| . The thick dashed lines 
show our adopted fitting formula. The figure shows that our 
fitting formula provides a good fit to the 'bright' sample, but 
significantly underpredicts the number of large EW systems 
in the 'faint' sample (Fig IA3|) . This is because our fitting 
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Figure A2. The observed Lyc* rest frame equivalent width 
(REW) distribution of the brightest 399 z = 3 LBGs of the sam- 
ple of Shapley et al. (2003) is shown as the black histogram. This 
sample had a median Mtjv = —22.8, and is labeled the 'bright' 
sample. The dashed line shows our adopted fitting formula (Eq|3j. 
Our fitting formula provides a decent fit to the data. 
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Figure A3. Same as Figure lA2l but for the faintest 398 LBGs. 
Our fitting function now underpredicts the number of large EW 
systems significantly. This is because our fitting formula adopts 
REW C = 22 A, while the best fit REW C = 33 A. Our fitting 
formula underpredicts the number of large EW systems. Applying 
a correction for this would lower our required values for F (see 
text), strengthening the conclusions of this work. 



formula adopts REW C = 22 A, while the best fit to this sub- 
sample is REWc = 33 A ( see the inset). On the other hand 
this choice is required in order to be consistent with data 
at fainter UV magnitudes (see below). We stress that our 
fitting formula underpredicts the number of large EW sys- 
tems at faint UV magnitudes. Applying the correction would 
lower our required values for F (see text), strengthening the 
conclusions of this work. 

• We choose the normalization^ factor such that all drop- 
out galaxies have — oi ^REW^REWmax- While there are 
drop-out galaxies with smaller REW, this choice automati- 
cally results in a fraction of drop-out galaxies with Lya in 
emission (i.e. REW > 0) that increases with Mtjv, as ob- 



This factor is g 



iven by N = [exp ( REW ° ( 1 il%) ) 



cxp 



-REW„ 



REW C (M UV ) j\ 



)1/REW C (M UV ). 
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Figure A4. This Figure shows the fraction of drop-out galaxies 
with a Lya REW > 50 A, as a function of Mtjv- The red circles 
represent the z = 3 data of Shapley et al. (2003, also shown in 
Fig [All and Fig EH). The blue squares represent data from Stark 
et al. (2010) from a Keck spectroscopic survey of z = 3 — 7 drop- 
out galaxies. The black solid line shows our fitting function. The 
black dotted lines show our fitting function when we increase or 



decrease the parameter 



by 4 A. Our fitting function is 



clearly consistent with the data (see text). 



served in the sample of Shapley et al (2003). This choice 
for the normalization constraint allows us to describe the 
Muv-dependence of the observed shape and normalization 
of the Lya REW PDF at REW> with one single parame- 
ter, namely REW C (we keep ai constant within this range). 
However, as we show next this single parameter description 
breaks down at fainter (Mtjv ;> — 21.5) magnitudes. 

• Stark et al. (2010) found that the REW-PDF is sensi- 
tive to Mtjv- This dependance is relatively weak in the range 
-22.0 < Mtjv < -20.5, but strong from -20.5 < Mtjv < 
— 18.5. This is shown in Fig IA4I where we compare the ob- 
served fraction of drop-out galaxies with REW> 50 A as a 
function of Mtjv to our fit. The red squares show the data 
presented by Shapley et al. (2003), while the blue circles 
show the data presented by Stark et al. (2010). The solid 
curve represents our adopted fitting function. We match the 
rapid evolution at Afuv > —20.5 by decreasing a\ as de- 
scribed in the paper. The dotted curves show the model when 
we increase or decrease the parameter dR ( ^j v by 4 A. Note 
that the data point at Mtjv = —21.9 lies a few a above 
our fitting function, because this data point lies significantly 
above the data at —21.5 ^ Afuv ^ —21.0. Forcing better 
agreement with this data point therefore automatically re- 
sults in worse agreement with the data at fainter UV mag- 
nitudes. This is why our best-fit model does not provide the 
best fit to the data shown in Figure [ 



A2 Alternative Parametrization of the EW-PDF 

To make sure our results do not depend on our chosen func- 
tional form, we also study an alternative parametrization 
of the Lya REW PDF. In this model, parametrize the ob- 
served shape and normalization of the Lya REW PDF at 
REW> at all Afuv with a single parameter (REW C ). 
To capture the observed evolution of the REW-PDF with 
Mtjv (weak in the range —22.0 < Mtjv < -20.5, and 
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Figure A5. Same as Figure IA4I but for our alternative 
parametrization of the REW-PDF (see text). 



strong from —20.5 < Mtjv < —18.5), we add a cubic term 
(AM) 3 into the expression for REW c (Mtjv)- Specifically, 
we keep the parameter ai = 20 A fixed at all Afuv, and 
REW c (Muv) = REW c ,o + i.(AM+[AM] 3 ) + ^Az. The 
data is well described by REW C , = 23±2 A and b = 4.4±lA 
(see Figure |A"5)) . 

We have repeated our analysis by fitting the observed 
Lya luminosity function for this alternative model. In our 
standard model (see § 14. 2|) we keep the evolution of the 
REW-PDF constant at Muv > -19.0. The results of this 
analysis are presented in Figures |A6I which shows that our 
constraints on F and dR ^ z Wc are very similar to the results 
already presented in the paper. We also found good fits for 
the UV luminosity functions. This gives us confidence that 
our results are not sensitive to the precise choice of the func- 
tional form that was used to model the REW-PDF. 



APPENDIX B: CALCULATION OF EW-PDF 
FOR LYa SELECTED GALAXIES 

The EW-PDF of galaxies as a function of Lya luminosity, is 



P(REW|L D 



dM uv P(REW|L Q , M uv )P(Muv|io 



Bayes theorem states that 



(Bl) 



P(M uv \L a ) = 



P(M uv ,L a ) P(L a \M uv )P(M v 



P(L a ) 



P{L a ) 



(B2) 



and we may write 
P(REW|L Q ) = 



dMuv P(REW\L a , M m )P(L a \Mvv) ^BTTT -( B3 ) 

We also know that P(L a \M vv ) = P(REW|M UV )^g^ = 
P(REW|M U v) s ^ (see § [2X2]) . For a particular combi- 
nation of L a and Afuv, the REW is fixed, implying that 
P(REW|L Q ,Muv) = 8o(g(M vv ), where 5 D (x) denotes the 
Dirac delta function and g(Muv) = REW — C xb.ewxl uv ■ 
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Figure A6. Same as Figure [T] but for our alternative parametrization of the REW-PDF. The most likely value for F decreased further 
from F = 0.43 to F = 0.34, and F = 1 is ruled out at >99% CL. 



We can therefore write 



/• — CD 

P(REW|L Q ) = / dM uv 8 (g(M uv ))P(M uv ) 
Jo 



REW 

;P(REW|M UV )- 



P(L 

a. 

P(REW|M uv ,c) REW 



ff'(Muv,c) L a P{L a ) 



P(M UV , C ) 



where in the last step we evaluated (/(Mtjv) at the pole of 
g(Mw), i.e. when C x REW x i(n - 4 ( M uv,c+s) = La We 
can simplify this further to 

P(REW|L Q ) = -2.5P(REW|M uv , c ) in f^; ( ^ ) . (B5) 

We finally compute the REW-PDF as 

P(REW) = Af J ° max dL a P(KEW\L a )P(L a ) 

= M [ La ^ P(REW|A/ U v,c)P(M uv ,c)dlog 10 i a (B6) 



where M is the normalization constant (which absorbed all 
numerical factors). 



© 2009 RAS, MNRAS Q00.[THT4l 



